Metal-insulator transition in the Edwards model 
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Abstract. To understand how charge transport is affected by a background medium and vice 
versa we study a two-channel transport model which captures this interplay via a novel, effective 
fermion-boson coupling. By means of (dynamical) DMRG we prove that this model exhibits 
a metal-insulator transition at half-filling, where the metal typifies a repulsive Luttinger liquid 
and the insulator constitutes a charge density wave. The quantum phase transition point is 
determined consistently from the calculated photoemission spectra, the scaling of the Luttinger 
liquid exponent, the charge excitation gap, and the entanglement entropy. 



The way a system evolves from a metallic to an insulating state is one of the most fundamental 
problems in solid state theory. Electron-electron and electron-phonon interactions are the driving 
forces behind metal-insulator transitions (MITs) in the majority of cases. For example, the 
Mott-Hubbard MIT pQ is caused by strong Coulomb correlations, whereas the Peierls MIT [2] 
is triggered by the coupling to vibrational excitations of the crystal. Theoretically the MIT 
problem can be addressed by the investigation of generic Hamiltonians for interacting electrons 
and phonons such as Hubbard or Holstein models [3]. In one dimension (ID), these models 
exhibit a MIT at half-filling, where on the insulating side of the MIT a spin-density-wave 
(SDW) or a charge-density-wave (CDW) broken-symmetry ground state appears, respectively. 
On the metallic side, near the MIT, charge transport then takes place within a strongly 
correlated "background" that anticipates the developing SDW, respectively CDW, order. Since 
the particles responsible for charge transport and the background order phenomena are the 
same, the problem is very complex. 

A path forward might be the construction of simplified transport models, which capture 
the basic mechanisms of quantum transport in a background medium in a rather effective way. 
Along this line a novel quantum transport model has been proposed recently [3], 
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This so-called Edwards model mimics the correlations inherent to a spinfull fermionic many- 
particle systems by a boson affected hopping of spinless particles (see Fig. [1]) . For the half-filled 
band case, the model describes a repulsive Tomonaga-Luttinger liquid (TLL), provided the 
excitations of the background are energetically inexpensive (ojq < wo,c) or W A1 readily relax 
(A > A c (wo))- This defines the fluctuation dominated regime. By contrast, strong background 
correlations, which will develop for large ujq and small A <C i& tend to immobilize the charge 
carriers and may even drive a MIT by establishing CDW long-range order [5]. 
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Figure 1. The Edwards model ([I]) describes a very general situation: As a charge carrier (•) 
moves along a ID transport path it creates an excitation with energy ujq (*) in the background 
at the site it leaves or annihilates an existing excitation at the site it enters. The background 
medium may represent, e.g., a magnetically, orbitally or charge ordered lattice. One assumes 
that the (de) excitation of the background can be parameterized as a bosonic degree of freedom. 
Any distortion of the background can heal by quantum fluctuations. Accordingly the A-term 
allows for spontaneous boson creation and annihilation processes. 



In the present work, we employ density-matrix renormalization group (DMRG) and 
dynamical DMRG methods [6] to analyse the ground-state properties of the Edwards model 
and the charge carrier dynamics for the limiting case of high-energy background fluctuations. 

Let us start with the discussion of the photoemission (PE) spectra. The single-particle 
spectral function probed by angle-resolved [inverse] PE reads 
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Here A~(k,u) [A + (k,uj)] is associated with the emission [injection] of an electron with wave 
vector k, i.e. fr = f k and f£ = fj,. \ipo) is the ground state of a iV-site system in the Nf- 
particle sector, while \ip^) denote the n-th excited states in the Nj ± 1-particle sectors with 
excitation energies = — Eq. For the half-filled Edwards model we have Nf = N/2. 

Figure [2] shows A(k, u) for a stiff background, i.e. the distortions induced by particle hopping 
are energetically costly. In this regime the bosons will strongly affect particle transport: The 
quasiparticle mass is sizeably enhanced and a renormalized band structure appears but — if A is 
large enough — the system remains metallic, as can be seen from the finite spectral weight at the 
Fermi energy Ep (left panel). As the system's ability for relaxation decreases, i.e., at fixed ojq, 
A falls below a certain critical value, a gap opens in the single-particle spectrum at &f = tt/2 
(middle panel). Evidently the system has become an insulator. We note the internal feedback 
mechanism: The collective boson excitations originate from the motion of the charge carriers 
and have to persist long enough to finally inhibit particle transport, thereby completely changing 
the nature of the many-particle ground state. The collective boson-particle dynamics leads to 
an asymmetric band structure for k < kp and k > kp (see inset). While the induced hole 
probed by PE can only move coherently by a six-step process with three bosons first excited and 
afterwards consumed, an additional electron can easily move by a two-step process even if strong 
CDW correlations exist in the background [5]. We note that the (I)PE spectra exhibit weak 




Figure 2. Line-shape of the single-particle spectral A(k,u)) in the half-filled band sector of the 
ID Edwards model. The insets shows the dispersion of the absorption/emission maximum. For 
the numerics, we consider an N = 32-site chain with open boundary conditions (BC) and map a 
boson site, containing 2 Ub states, to rib pseudosites. We take up to 4 pseudosites, keep m = 500 
density-matrix eigenstates, and use a broadening r/ = 0.1. All energies are given in units of tb- 



signals around the bare boson energies ±cjo (not shown). The interrelation of charge dynamics 
and background fluctuation becomes apparent again, if we decrease uq keeping A fixed (right 
panel). Now the fluctuations overcome the correlations and the system returns to a metallic 
state which is different in nature, however, from the state we started with: A(k, uj) shows sharp 
absorption features near kp only and is "overdamped" at the Brillouin zone boundaries, where 
the spectrum is dominated by bosonic excitations. 

In order to determine more precisely the phase boundary between the metallic and insulating 
ground states, typifying a Tomonaga-Luttinger liquid (TLL) and a CDW, respectively, we 
analyse the limiting (N — > oo) behaviour of the TLL charge exponent 

K p = nKm^, with 5c (g) = i^e i ^- fc )((n i -i)(n fc -i)), q = ^, (3) 

as well as those of the single-particle (charge) gap, A C (N) = Eq + Eq — 2Eq, and monitor the 
finite-size scaling of the entanglement entropy difference [7] 



AS N = S N (N/2) - S N (N/2 - 1) = -— In cos — , (4) 



3 



7T 



N 



where Sn{1) = — Tr[p; Inpi] = ^- In 
decreases from K p = 1, as A is lowered, and fina 
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+ s\. We expect that the TLL charge exponent 
ly reaches 1/2 at the MIT point, if the transition 



is of Kosterlitz-Thouless type [H [9] . The central charge c* should scale to unity in the metallic 
TLL regime [TO] . 

Figure [3] demonstrates that the iV — > oo extrapolated K p indeed becomes 1/2 at some critical 
value, where A~ 1 (wo = 10) — 5.89, indicating the MIT. In the metallic phase we find a repulsive 
particle interaction, K p < 1. Our DMRG results point towards an exponential opening of the 
charge gap entering the insulating state, which corroborates the Kosterlitz-Thouless transition 
scenario. Note that the CDW state of the Edwards model is a few boson state, in contrast to 
the Peierls CDW phase of the Holstein model [5] . That means the MIT in the Edwards model is 
driven by strong correlations, as for the Mott-Hubbard transition. To extract the central charge 
c* we use the entanglement entropy difference, Eq. (U|), rather than directly exploiting Sjy(l). 



Figure 3. Left panel: N —¥ oo extrapolated value of the TLL parameter K p , respectively of 
the charge gap A c , as a function of A -1 for ujo = 10 (open BC). Middle panel: Entanglement 
entropy differences ASV for different system sizes (periodic BC). The red dashed-dotted line 
gives the MIT transition point in reasonable agreement with the value obtained from K p (black 
dashed line). Right panels: Critical value of A" 1 (filled circle, top panel) and central charge 
c* ~ 1 (filled square, bottom panel), both extrapolated from the maxima of ASjv- Here we use 
m = 2000, rib = 2, and ensure a discarded weight less than 10" 10 . 



For a model with spinless fermions this is advantageous because we can work with a fixed system 
size, thereby avoiding antiperiodic BC that give rise to complex phase factors [7]. As can be 
seen from the middle panel of Fig.O for A" 1 < A^ 1 , the rescaled quantity — 3ASn/ ln[cos(7r/iV)] 
extrapolates to unity as N — > oo. This opens an alternative route to detect the MIT point. We 
find that the A~ 1 (wo) determined by extrapolating the maximum of ASn, i.e. in a completely 
different manner, matches the critical value obtained from K p surprisingly well. Simultaneously, 
indeed c* — >• 1 (see right panel). 

To summarise, we have studied the spectral and ground-state properties of the ID Edwards 
fermion-boson transport model by large-scale (dynamical) DMRG numerics. We showed that 
strong correlations within the background medium will not only affect the charge-carrier's 
dynamics by enhancing the quasiparticle mass but may even trigger a metal-insulator quantum 
phase transition. The MIT transition point has been determined in good agreement both from 
the TLL charge exponent and the entanglement entropy difference. We stress that to date only 
a very small number of microscopic model exists which have been rigorously shown to exhibit a 
MIT. 
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